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Abstract 

In this paper we study the fluctuation spectrum of a linearized 
Vlasov-Poisson equation in the presence of a small external electric 
field. Conditions for the control of the linear fluctuations by an exter- 
nal electric field are established. 



1 Introduction 

In the past, the fluctuation spectrum of charged fluids was studied either by 
the BBGKY hierarchy derived from the Liouville or Klimontovich equations, 
with some sort of closure approximation, or by direct approximations to the 
N-body partition function or by models of dressed test particles, etc. (see 
reviews in p] [2]). 

Alternatively, by linearizing the Vlasov equation about a stable solution 
and diagonalizing the Hamiltonian, a method has been developed [3] that 
uses the eigenvalues associated to the continuous spectrum and a canonical 
partition function to compute correlation functions. Here this approach will 
also be followed to study the control of the fluctuations. For simplicity we 
will consider the one-space dimensional case. 
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A Vlasov-Poisson equation with an external electrical field control term 
is considered. Following the method developed by Morrison [I] we use an 
integral transform to solve the linearized equation. With a view to applica- 
tions to more general kinetic equations (gyrokinetic, etc.) we also discuss in 
the appendix a generalization of Morrison's integral transform. 

Control of the Vlasov-Poisson equation 

9tf + v- Vx/ + Vx(p ■ Vvf = C (t, x, v) 
in a periodic domain 

(t,x,v) G [0,T] xT n xR n 

by means of an interior control located in a spatial subdomain has been 
discussed by Glass j5]. Conditions for controllability between two small dis- 
tribution distributions fo and f\ were established. However, to steer f to fi 
a control C (t, x, v) that depends on the velocities is required and it is not 
clear how such a control could be implemented in practice. Therefore, we 
have restricted ourselves to the more realistic situation of a (small) control- 
ling external electric field. In addition we concentrate on the problem of the 
damping of the small oscillations around an equilibrium distribution. 

In Sect. 2 the linearized Vlasov-Poisson equation with control is solved 
by an integral transform and in Sect. 3 two controlling problems are studied, 
namely the control of the total energy of the fluctuations by a constant 
electric field and the dynamical damping of the fluctuating modes by a time- 
dependent electric field. 



2 The linearized equation with control 

Consider a Vlasov-Poisson system in 1 + 1 dimensions 
df , _df e {d<$> f „ cl __ ^ df 



dt V dx m\dx ^ ( X ^)J q v ® (-0 

—® f = -e J f( v )dv + p B (2) 

with an external (control) electric field E c (x, t) and a background charge 
density ps (x) chosen in a such a way that the total charge vanishes. Consider 
now the linearization about a homogeneous solution. 

f(x,v,t) = f<®(v)+6f(x,v,t) (3) 
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with 

such that 
Then 



A$f = ~ e J f m ( v ) dv + Pb 



A*f = = V *f 



and (v) is indeed a homogeneous static equilibrium. The linearized equa- 
tion is 

05 f 05 f edSQfdfW e ^. , J/< ' 
at ox m ox ov m ov 
where we have assumed that 

E c (x,t) = 0(5f) (5) 

that is, E c (x,t) is a small external (control) electric fielcfl. 
Fourier transforming all the perturbations 

5f(x,v,t) = Y,fk(v,t)e ikx (6) 

5$ (X, t) = fa (*) ^ (7) 

E c (x,t) = ^i? fc (i)e^ (8) 

leads to 

d t f k (v, t) + ikvf k (v, t)-i^-(f f k (ji) dfx) /<°>' + -£ fe (t) /(°)' = (9) 

With a view to applications to more general kinetic equations (gyrokinetic, 
etc.) a more general equation is studied in the appendix. Equation (jUJ) is 
then a particular case of Eq. fl24|) with 



g 1 (v) = ikv, gt ( v ) = -i*-fW(v), gs(y) = l, C (v,t) = -E k (t) f®' (v) 

(10) 



1 that is, the control electric field is of the same order as the fluctuations, not of order 
J Sf (v) dv, which would lead to a trivial control situation 
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the integral transform being (as in Morrison [1]) 

G k (it, t) = (G.A) («) = (l - ^ (V 0) ')) /* ^+^f m ' W*) (11) 
with left inverse 

°--{( 1 -£(*'*)H£'" p )T (12) 

G + being 

(G+A) («) = ( i - ^ (#/ (0) ')) /* («) - ^/ (0) ' W*) (is) 



Applying the integral transform (ITTj) to Eq.Q it becomes 

<9 t (GL/) («) + zfcu (u) = (t) / (0) ' («) (14) 

with solution 

G fc (u, t) = e"^ f G fc (u, 0) - ^/W' (it) jf £ fe (r) e^dr) (15) 
Then, according to (|T2|) . the Fourier modes solution is 

/fc(M) = — 4 2 ( ie ) 

(i-S (*)) +(Sf {oy (v)) 2 

™ 2 f (0) ' M# f Gk ( "' t] I (r /) 

(1 („)))» +(^/(W („))»' 



m/c 2 



2.1 Control of the linear modes by the electric field 

Nonlinear stability of the steady states of the Vlasov-Poisson equation when 
the phase-space density is a decreasing function of the particle energy or 
depend on other invariants has been studied[6] [7] by the energy- Casimir 
met hod [8]. This means that deviations from the steady-state will remain 
bounded in time. 
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However, as expected from the non-dissipative nature of the Vlasov equa- 
tion, the linear fluctuation modes of the uncontrolled equation are oscillatory 
and, once excited by a perturbation, they will not decay. As shown by Morri- 
son m they may be used to obtain a statistical description of the fluctuations 
by the construction of a partition function. Here, one focus on the control 
of the fluctuations by the external electric field. Two situations will be con- 
sidered. The first considers a constant in time electric field and tries to 
minimize the total energy associated to the fluctuations. The functional to 
be minimized is 

fT 

F 1 (E k )=\im / dtdu\G k {u,t)\ 2 (17) 

In the second situation we allow the electric field to be time-dependent and 
chosen in a such a way as to introduce a damping effect in the solution (|T5|) . 

For the first case (E k independent of time), with the solution (fT5|) one 
obtains a minimum for the functional Fi at 

- / g|g (-ImGfc (u, 0) + iReG k (u, 0)) du 
^-f(^) 2 du 

mk J \ u I 

For this electrical field F\ is 

i 71 imin = / ((ReGfc (u,0)) 2 + (lmG k (u,0)) 2 ) du 

J f f(°)' (u)ReG k (u,0) \ + f f^' (u)ImG k (u,0) \ du j (18) 
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2/1 ^M] d u 



a smaller value as compared to the case E k = 0, which would be Fi (E k = 0) = 
/ ((ReG k (u, 0)) 2 + (ImG fc (u, 0)) 2 ) du. 

In the second case one allows the electric field to be time-dependent. 
One aims at controlling the fluctuation modes by an electric field induced 
dynamical damping. One looks for the solution of 



lim (g h (u, 0) - -/<°>' (it) / E k (r) e iTku dr) = (19) 



obtaining 



mk r G k (u,0)^ lku , 



E k (t) = — / K \ ' > e~ lkut du (20) 
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Then with this electric field 

-t 



G k (u, 0) - -/W' (it) / E k (r) e iTku dr (21) 



G k (u,0)-—f^ (u) / du 



^ikyu—u jt 



2vr J v V-oo / (0) V) ifc(u-u') 
and from 

, ik[u—u )t , 

h p \ J — 1 / A 

5[u-u) (22) 



27T ik (u — u) t^co 

one sees that the electric field ( J20l) induces a dynamical damping of the 
fluctuation modes. 



3 Appendix. An integral transform for lin- 
earized kinetic equations 

Morrison [I] solves the linearized Vlasov-Poisson equation by a Hilbert trans- 
form. However, for some practical applications, the linearized kinetic equa- 
tions are more complex. For example the gyrokinetic Vlasov equation written 
in gyrocenter phase-space coordinates is [9] 

f + i. v/+( 7§ = (23) 

where U 

This is the motivation to study an equation more general than the lin- 
earized Vlasov-Poisson (jUJ). Linearized Fourier kinetic equations are of the 
type 

^l + gi (v)f(v)+g 2 (v) JgaQi)f(li)dn + C(v,t) = (24) 
with gi a monotone function of v. Let T be a transform such that 



(TgJ) (u) = J g 3 (p) / (//) dp + g x (//) (Tf ) (u) (25) 

namely 

(Tf)(u) = p[ 9 ^ )nv) ( dv (26) 
J gi(v)- gi (u) 
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Notice that, for invertible gi, the T— transform may be written in terms of 
the Hilbert transform 

(Tf) (u) = tt (h 9 4- o gT 1 W («)) (27) 

Then, one defines 

(«) = (1 + (Tft) («)) / («) + <? 2 («) (T/) («) (28) 

and 

(GL/) (tt) = (1 + (Ta.) (u)) / («) - <? 2 («) (T/) (tt) (29) 
|(1 + Tg 2 f + vr 2 (^) 2 | G_ is a left inverse of G + 

(G-G+f) (tt) = <!(!+ T^ 2 ) 2 + vr 2 (^p) 1 / (tt) (30) 



as may be checked using (127]) and the properties of the Hilbert transform 

Din- 

(^--transforming Eq. (|24|) one obtains 

9 t (G_/)+ 5l («)(G_/) = -C(u,t)+g 2 (u)T(C)(u)-C(u,t)T(g 2 )m) 

= l(u,t) 

with solution 

G_ (tt, i) = e-* Sl(u) (g- (tt, 0) + jf 7 (tt, r) e T5l(u W^ (32) 
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